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We describe the fully gauge invariant cosmological perturbation equations in teleparallel 
gravity by using the gauge covariant version of the Stewart lemma for obtaining the vari¬ 
ations in tetrad perturbations. In teleparallel theory, perturbations are the result of small 
fluctuations in the tetrad field. The tetrad transforms as a vector in both its holonomic 
and anholonomic indices. As a result, in the gauge invariant formalism, physical degrees 
of freedom are those combinations of perturbation parameters which remain invariant 
under a diffeomorphism in the coordinate frame, followed by an arbitrary rotation of the 
local inertial (Lorentz) frame. We derive these gauge invariant perturbation potentials 
for scalar perturbations and present the gauge invariant field equations governing their 
evolution. 
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1. Introduction 

Standard cosmology, based on the assumption of isotropy and homogeneity at large 
scales, has remarkable success in describing the known properties of the universe. 


1 
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The origin of local inhomogeneous structures such as galaxies and galaxy clusters 
then can be traced to small fluctuations around the FRW background of standard 
cosmology. Generally in cosmology, computing the perturbation of a quantity in¬ 
volves determining the difference between the value of this quantity in the actual 
physical spacetime and the value of it in the reference unperturbed FRW frame. 
To determine this difference, it is necessary to compare the values at the same 
spacetime point in both frames. To do this, one needs to choose a map that shows 
which point on the perturbed spacetime is the same as a given point in the back¬ 
ground geometry. Choosing a map is called a gauge choice and therefore changing 
the map implies performing a gauge transformation. Generally the choice of the 
map is arbitrary. This is referred to as the gauge freedom of the perturbation the¬ 
ory. There are two different approaches in performing calculations in this setup. 
One method is to choose a gauge (or identification map) and perform the calcula¬ 
tion in this gauge. This method, while it’s usually simpler, has a drawback. It may 
result in the so called ’gauge mode’ solutions which are not real physical modes. 
The second method, first introduced by Bardeen in) 1 involves finding some gauge 
invariant combinations of perturbation modes and rewriting the equations in terms 
of this gauge invariant quantities. This approach has the advantage of having only 
real unambiguous physical quantities but it is technically more involved. The gauge 
invariant theory of cosmological perturbations in general relativity, has been devel¬ 
oped through the years and applied to various problems including inflation theory 
and computing the spectrum of the cosmic microwave background radiation. For 
thorough reviews of this method see references 

Teleparallel theory of gravity first introduced by Einstein® in an attempt to 
unify gravity with electromagnetism. This theory can be regarded as the gauge 
theory for the translation group and in its general form does not possess the same 
symmetrical structure as the theory of general relativity (GR). 7 Translational gauge 
theory is formulated in terms of coframe: in each point of a n-dimensional manifold 
one introduces n linearly independent vector as the basis of spacetime and their dual 
covectors or coframes. The theory then possess a torsion field which can be regarded 
as the translational field strength corresponding to the coframe field P The coframe 
theory in its general form, does only possess the diffeomorphism invariance and the 
invariance under the global Lorentz transformation of tetrads. It is not invariant 
under a local Lorentz transformation. Restoring the local Lorentz invariance means 
imposing some restrictions on the form of the Lagrangian. Doing this will result in a 
theory which is dynamically equivalent to general relativity and is usually called the 
teleparallel equivalent of general relativity (TEGR) in the literature.^ 10 Teleparallel 
gravity and its extensions 11 12 have generated renewed interest in recent years, in 
the hope that it may offer solutions to some of the cosmological problems like the 
origin and the nature of the dark energy P-^31 It should be stated here that both 
teleparallel gravity and general relativity can be regarded as special cases of a more 
fundamental gauge theory of gravity called Poincare gauge theory (PGT) which 
contains both torsion and curvature as translational and rotational field strengths 
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respectively . 15 Teleparallel gravity could not achieve the unification of forces desired 
by Einstein, nonetheless there are significant incentives in studying this theory and 
its properties. For example, according to some authors a quantization approach 
based on the teleparallel variables will probably appear much more natural and 
consistent compared to general relativity P Also it has been shown that owing to 
its gauge structure, attempts to incorporate the nonlocal effects into the theory of 
gravity can best be done in a teleparallel formulation . 16 

For the reasons stated above, it seems interesting to study the cosmological 
consequences of teleparallel theories of gravity. One of the differences between GR 
and TEGR is in defining the dynamical variables of the theory. In TEGR unlike GR 
tetrad fields play the important role of dynamical variables. This fact has motivated 
the authors to consider the cosmological perturbations by directly perturbing the 
tetrad fields and study the potential differences. In this paper we also attempt 
to obtain fully gauge invariant formulation of equations governing perturbation 
dynamics. 


2. Notation and definitions 

Throughout the paper, the Greek indices ■■ run over 0,1, 2, 3 and refer to the 
spacetime coordinates; The Greek letters from the beginning of the alphabet a, /?,... 
run over 1, 2 , 3 and refer to spatial coordinates, middle Latin letters i,j ,... run over 
0,1, 2,3 and refer to the 4D tangent space coordinates and finally Latin letters from 
the beginning of the alphabet a,b ,... run over 1, 2, 3 and indicate the spatial tangent 
space coordinates. 

As stated, in teleparallel gravity one considers a set of n linearly independent 
vectors which form a basis in the tangent space on every point of the 

manifold. The dual of this basis = e l ,/lx 11 are coframes. The dynamical variable in 
TEGR (e z M ) are called the tetrads and they relate anholonomic tangent space indices 
to the coordinate ones. The spacetime metric is not an independent dynamical 
variable here and is related to the tetrad through the relations 


(1) 

■CO 3. 

QJ 

II 

a. 

(i) 

— Tjij 

(2) 


The inverse of the tetrad is defined by the relation ef = 6?. 

Teleparallel geometry, X 4 is obtained by the requirement of vanishing curvature. 
In the special case of TEGR the spin connection of the theory is also assumed to 
be zero. This assumption is usually called the absolute parallelism condition and 
by imposing it the connection of the theory will be Weitzenbock connection defined 

a:P 


pp _ p po J 
1 fiv ~ U v^n 


( 3 ) 
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which unlike Livi-civita connection is not symmetric on its second and third indices. 
The curvature of this connection is identically zero and the torsion tensor is 

T% v = e?{d,ei - d ^). (4) 

Contorsion tensor which denotes the difference between Livi-civita and Weitzenbock 
connections is 

K^ p = -\{T^ p -T^ p -T p n (5) 

and the superpotential tensor is defined as 

Sr = + %T av a - a) • ( 6 ) 

In correspondence with Ricci scalar, one can define torsion scalar 

T = S p f T p f (7) 

The gravitational action in TEGR is 

/= i6b/ ± ' |e|r (8) 

where |e| is the determinant of e“ and from the relation (1), one can easily finds it 
equal to \]~g. Variation of the above action with respect to tetrads gives the field 
equations of TEGR 

e ~ 1 d l _i{ee f ’S p n - e*T p pX S p ^ + \e\T = 4ttG efs; (9) 

where H. 1/ is the energy-momentum tensor. 


3. Gauge invariant cosmological perturbations in teleparallel 
gravity 

We now begin the process of deriving fully gauge invariant cosmological perturbation 
equations in teleparallel gravity. For a given quantity Q , its perturbation 6Q is the 
difference between its value in the actual perturbed spacetime and in the background 
reference geometry, calculated at the same point 

SQ(p) = Q(p)-Q°(D- 1 (p )) (10) 

where D is the diffeomorphism D : N —> M from the background manifold N to the 
physical manifold M. Another diffeomorphism, say D , result in a different SQ(p) 

6Q(p) = Qlp)-Q°(D- 1 (p)) (11) 

The change of the identification map between the two manifolds D — » D results 
in a gauge transformation 6Q(p) 5Q(p). As another viewpoint, a gauge change 

can also be regarded as a coordinate transformation on the background manifold 

s'* -)• = a? + ^ (12) 



May 31, 2016 0:34 WSPC/INSTRUCTION FILE GICP 


5 


This transformation will lead to a change in SQ as 

A 8Q = SQ-SQ = L^Q (13) 

where L is the Lie derivative in the direction of £. A quantity then called gauge 
invariant if 


A SQ = Li-Q = 0 (14) 

This important result was first derived irPO (see alscPD) and is usually called the 
’Stewart Lemma’ in the literature. 

As mentioned before teleparallel theory in its general form can be regarded as 
the gauge theory for the translation group. In such a gauge theory, the notion of 
ordinary Lie derivative used in Stewart lemma in general relativity, will no longer 
be adequate. Here, the variation of a field is obtained by using the gauge covariant 
Lie derivative with respect to a vector v, defined for a Lie-algebra valued form T 
apMH 


L®* = V \D^> + D(v\ V) (15) 

where J is the interior product and D is the covariant derivative, here obtained 
by using the Weitzenbock connection (7). In coordinate notation, using the gauge 
covariant form of the stewart lemma, the variation in the perturbations of the tetrad 
field will be 


A<5e* 




stv„r 


(16) 


4. Tetrad perturbations and gauge invariant potentials 

We begin by considering the unperturbed Friedmann-Robertson-Walker line ele¬ 
ment using a conformal time parameter 

ds 2 = a 2 (r){— dr 2 + y a p dx a x^} (17) 


where y a p is the metric of the spacelike surfaces of constant torsion in teleparallel 
gravity and we denotes the covariant derivative associated with it by D a . This metric 
will be used to raise and lower indices in the spatial hypersurface, for simplicity we 
assume yP = SP where fiP is the Kronecker delta. 

The simplest tetrad describing this FRW background is given by 


(a 0 0 0\ 
0 a 0 0 
0 0 a 0 
\0 0 0 a) 


As we are working in a TEGR setup, any other possible FRW tetrad is related to 
this by a Lorentz transformation and will result in the same dynamics. 
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The tetrad transforms as a vector in the coordinate space so generally like any 
other vector A a , they can be decomposed to a scalar part and a purely vector part 
according to 

A a = D a (j> + B a (18) 

where 0 is a scalar field and B a is a solenoidal vector i.e. D a B a = 0. The general 
perturbed FRW tetrad of teleparallel gravity is given by 

e 0 = “k 4l) > &a = a Cot ) Gq = a B a , 

e a a =a(6 a a +D a a ) (19) 

and its inverse 

hl = -{l-A) , b% = —B° , h° a = --C it 

a a a 

K = l (<£ - D“) (20) 

where an over-dot denotes a tangent space index. Here A is a scalar and B l , C a 
and D l a act as vector degrees of freedom. The tensor part of the perturbations 
in this formalism, comes from the metric as can bee seen in the appendix where 
D af 3 = VijS^Dp contains the tensor part of the perturbations. Note that D a p is still 
of the first order in tetrad perturbations. 

If we decompose in (12) into temporal and spatial parts as 

= (T,L a ) (21) 

then using stewart lemma (14) and using the covariant form of the Lie derivative 
defined in (15), under transformation (12) we have 

A A = T' + hT , A B i = L h + hU (22) 

AGq, = 8 a T , A D l a = D a V 

here a prime denotes differentiation with respect to proper time and h = a'/a. For 
scalar perturbations, we should have 

R ?; = ITS , C a = D a C , D“ = D (23) 

the scalar parameters transform as 

H A + T' + hT (24) 

B B + L' + hL 
C^C + T 
D^D + L 

where A = D a D a is the spatial Laplacian. The above perturbation parameters are 
not gauge invariant; however there exist two different combination of them which 
are gauge invariant 


= A - C' - hC 
$ 2 = B-D'-hD 


( 25 ) 
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These two parameters are gauge invariant potentials and corresponds to Bardeen’s 
potentials of general relativity. It is obvious that they are not equivalent. By rewrit¬ 
ing the perturbed field equations in terms of these parameters, one can deal with 
only real physical quantities and any gauge ambiguities will be removed. The reason 
for different gauge invariant potential in TEGR and GR is twofold. One reason is 
that in teleparallel gravity, one should use the teleparallel version of the Lie deriva¬ 
tive in the Stewart lemma (14), defined in a weitzenbock T 4 geometry. The other 
source of difference is the fact that unlike metric in general relativity, tetrads does 
not have any tensor modes. 

5. Matter perturbations 

We assume an unperturbed energy-momentum tensor in the form of a perfect fluid 
given by 

= (p + p)u^u v + pg^ (26) 

where u M is the 4-velocity. Under the gauge transformation we have 

S^ap —■t S^ a p + Tpa 2 pS a p + 2pa 2 D( a Lp) (27) 

fcoo —t fc 0 o T poi 2 [pT + 2 T' + 2hT\ 

S^oa —t <5^oa + a2 [pL' a + ph.L a + pD a T] 

where p = ^ and P = r ~- The perturbation in the energy-momentum tensor can be 
written as 

SE^ = pSpu fJ, u l/ + 2 aq^u v) + pSp(uV + O + pa 2 n^ (28) 

where is the momentum perturbation and is the anisotropic stress. Under 
gauge transformations we have 

Sp 5 p + pT + 2 T’ + 2 hT 

q a q a - pL' a - pD a T - ph,L a 

Sp —> Sp + pT + -D a L a 

o 

2 

TTct/3 ^ 7T a /3 T — SapD^L^ (29) 

For scalar perturbations L a = D a L, q a = D a q and Tr a p = A a pn , so 

Sp -> Sp + pT + 2 T’ + 2 hT 

q —> q — pL' — pT — ph.L 
2 

Sp —> Sp + pT + -A L 

o 

7T —^ 7T T 2L 


(30) 
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There are four different gauge invariant scalar combinations of matter and tetrad 
perturbation variables 


= Sp - 2A - pC 

&q = q + pC + pB 

3> p = Sp - pC - D 

= 7T - 2 D (31) 

In the next section we will present various components of the teleparallel field 
equation (9) in terms of the six gauge invariant variable presented here. 


6. Gauge invariant field equations 

The torsion, contorsion and superpotential tensor of the perturbed tetrad (19) is 
given in the appendix. Various components of the teleparallel field equation are 


(i = 0,i/ = 0) : 


1 






9 h 2 

d a B a ) + 2^ (1 “ 3A) =47rGe ' 


a. -^0 
O'— 'a. 


(32) 


(i = a, v = 0) : — 


2 a 3 


- d (a B a ) + S^dpB 0 + dl a C a] + D {a a) - 620% 


^6Z(A-\D\) + ^D\ 


a 

AirGe' 


—a n 0 --—r 

a 5 10 °] 2 a 5 a 


a 77 0 

a '—’a 


(33) 


(i = a,v = a) : dp 


- d [a D 0 a] + d [0 D a a] + d [0 D a “l + 25 a 0 d [a D v v] - 2 5 a a d [0 D r ' v] 


—a , 

~ 2 ~ 


( - d {a B^ + dZdpB 0 + - 8°D% + £>'<“ + dl°C a] 






d {a B a) 


a z 


^SZdpB 0 - 

h 
2 a 3 


a' 


£>' ( “ - -jD' 
or 


h ' 

r °) 

5 h 2 
2 a 3 


+ 


§£«(1 - - V3 al " c -] ~ 


(34) 
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(i = 0, v = a) 




£)'{&<*] _ — 


—B a + ad [a D 0 
l a 




a*°° 

3d^ a D 0 0] +d^D 0 a] 


d [0 C° 


\ha{C° 
5 h 2 


2 a 5 


B° 


-B°)\ 

= 4tt Gegs^ 


(35) 

These equations can be brought to the gauge invariant form with the use of the gauge 
invariant variables defined in (25) , (31) and the background field equations.After 
some manipulations, the gauge invariant field equations are 


1 / h' + <z. . 

_ 4 2 + _ (5 . 1 + 4 . 2 ) 


5?7 2 

5? + -^d v d^ 2 - $!) 


4?rGa I $*• - 2 h$ q + $ p 1 5% 


(36) 


f M>', - $' 2 - j = 4^G^($ P + $,)«*' 


—3h h h — 1 . 

o A4> 2 + ^t($ 2 - -2$l) = 4 ttG —<P p 
or or a z \ a 


(37) 

(38) 


^9 0 ($2 + /i$i) = 47rGa($ p + $ g )w a (39) 

where the first two equations are dynamical equations and the last two will act as 
constraints. Note that u v and Ui are the 4-velocity of the unperturbed background 
energy-momentum tensor. 


7. Dynamics of an Scalar Field as an example 

In order to find the observational consequences of the formalism described above, 
we examine a case where the matter content is described by a single scalar field. 
This can be applied to the inflationary era or dark energy models. In this case, for 
the background values of the energy - momentum tensor we have 

P=\j> 2 + V{4>) , P = (40) 

Perturbing the scalar field as 

<j)(x, t) = <j>(t) + S<j>(x, t) (41) 

and substituting in eq. (28), we get the components of the perturbed energy mo¬ 
mentum tensor. 
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As usual we expand an scalar fields x into its Fourier modes 

X(t,x)= f -^ 3 -\k{t)e lk S 
J (27r) 2 


(42) 



Fig. 1 . Scalar teleparallel perturbation potentials 4>i and its general relativistic coun¬ 
terparts for chaotic potential V (<j>) = ^ m 2 <f > 2 . The plots are for the wavenumber value 
k = yjjQQ hMpc~ l . The time axis is in Gigayears 


Here we are concerned with a case where there is no anisotropic stress and no 
momentum perturbation. In this case we can solve equations (39) and (40) together 
to find the two invariant perturbation potentials and $2 and then use one of the 
constraint equations (41) or (42) to find the gauge invariant density perturbation 
4>p. The numerical analysis is done for three types of cosmologically interesting 
potentials. 

• Chaotic potential: Figures (1) and (2) show the time evolution of geometric 
potentials $1 and $2 versus the cosmic time when the scalar field potential is 
assumed to be of the chaotic type i.e. V(4>) = \rn 2 (f > 2 . For comparison, we also 
plotted the corresponding potentials deriving from the theory of general relativity. 
As can be seen from the figures, the results are different at early times (high energies) 
but coincide at late times. 

• Exponential potential: Figures (3) and (4) show the time evolution of geometric 
potentials $1 and $2 versus the cosmic time for potential in the form of V{(j>) = 
Vbe _Al K Here again the results show the same trend, however the potentials coincide 
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Fig. 2. Scalar teleparallel perturbation potentials <I >2 and its general relativistic counter¬ 
parts for the chaotic potential = \rr? (f> 2 . The plots are for the wavenumber value 
k = yjjQQ hMpc~ l . The time axis is in Gigayears 


at at a slightly later time than the previous case. 

• Power law potential: Figures (5) and (6) show the results for potential in the 
form of V(4>) = A</> 4 . In this case the potentials for teleparallel and general relativity 
theories coincide at an earlier time than the first case. 



Fig. 3. Scalar teleparallel perturbation potentials 4>i and its general relativistic coun¬ 
terparts for exponential potential = Vge . The plots are for the same value of 

wavenumber k as in previous figures 
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Fig. 4. Scalar teleparallel perturbation potentials <f >2 and its general relativistic counter¬ 
parts for exponential potential V{(j> ) = Voe - *^ and the same value of the wavenumber k 
as before . 



Fig. 5. Scalar teleparallel perturbation potentials <f>i and its general relativistic counter¬ 
parts for potential V{(j)) = A0 4 . Other parameters are the same as previous figures. 

8. Conclusion 

The theory of TEGR is obtained from gauging the translational part of the Lorentz 
group and after applying some restrictions on the PGT Lagrangian will turn into 
a theory which is also invariant under local Lorentz transformation. These prop¬ 
erties has led the theory to be dynamically equivalent to GR while they still have 
some conceptual differences. This equivalence has led to a reluctance in studying 
many physical issues in the context of TEGR but the special geometric structure 
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Fig. 6. Scalar teleparallel perturbation potentials <f >2 and its general relativistic counter¬ 
parts for potential V(cp) = A<jf> 4 . 


of this theory has some interesting features to study. In the context of cosmological 
perturbations the results of directly perturbing FRW vierbeins, are found to be dif¬ 
ferent from the results obtained from perturbing FRW metric in general relativity. 
This difference stems from different transformational properties of tetrad and metric 
under a gauge transformation. These transformational properties are obtained by 
considering Lie derivative of the dynamical variables of the theory. In a teleparallel 
setup, one may use the gauge covariant Lie derivative which is not equivalent to 
the ordinary Lie derivative used in general relativity. Another difference comes from 
the fact that tetrad perturbations only have scalar and vector modes in the coor¬ 
dinate frame. As a result, the tensorial part will manifest itself as a combination of 
the vector perturbations modes in the metric. Studying cosmological issues in this 
context and comparing the results by observational data, concerning all the points 
mentioned above, seem to be interesting. In this paper also the gauge invariant 
perturbed field equations have been found and by means of gauge potentials. The 
formalism developed here can be readily extended for use in f(T) theories of grav¬ 
ity. However in the f(T, case extreme care is needed in defining the gauge invariant 
potentials as some of these models show violations of local Lorentz invariance. In 
that case a gauge invariant potential should remain unchanged under both a general 
coordinate transformation and a local lorentz rotation of the frame. 


Appendix A. 

In this appendix we present the non-zero components of the torsion and contorsion 
tensors, superpotential and the torsion scalar of the perturbed tetrad (19) 


Metric and inverse metric 
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900 (1 -|- 2A) , 90a ® (-^a ^ot) ? 9a(3 ^ (^a/3 “1“ -^a:/?) ? 

5 00 = -4(l-2^) , 5 °“ = 4(C a - sa ) , g aP = 2-(S° p - D aP ) 

a z a z a z 

Weitzenbock connection coefficients 

r\p = de D1 a , r7 0/ 3 = ^-B 7 , r^ 0 = M 2 + J D\ , r 7 00 = £' 7 

r° a/3 = ^c a , r ° 0/ 3 = dpA , r ° a0 = c' a , r \ 0 = h + A' 

Torsion 

T° a =C' a ^d a A , T°p = 2d [a C/3] 

T| 0 = W|-a^B“ + D^ , Tjfy = 29 [| g£)“ ] 

Contorsion 

RT“° = zi[c7'“ -d a A 

ar J 

4 |2£>' [a/3] - - d [a C 0] 

a z 

K a0 p = 4 [ - + d&Cfl + h( 1 - 2 A)6% - D a p 

1<a i = [ d[0c>a j] - d [a D\ } - d [a D 7 « 

+2h{6 1 [a C® - 5j a Bfl) 

Superpotential 

5 7 af> = 2- [ - c^ZT 7] + + 2<5“c^£’ 7 ??] 
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2S^d [a D v 7i] - 2(5j a C'ft - 6j a dflA) 
+ 2h{5j a C® - 5J a Bft) 



Torsion scalar 


T = -64 dnBP 

a 2 
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